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We consider a real massless scalar field in 1+1 dimensions satisfying time-dependent Robin bound- 
ary condition at a static mirror. This condition can simulate moving reflecting mirrors whose mo- 
tions are determined by the time-dependence of the Robin parameter. We show that particles can 
be created from vacuum, characterizing in this way a dynamical Casimir effect. 
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I. INTRODUCTION 



The phenomenon of particle creation from quantum 
vacuum by moving boundaries or due to time-dependent 
properties of materials, commonly referred to as the dy- 
namical Casimir effect (DCE) [1, 2], has been investi- 
gated since the pioneering works of Moore [3] and DeWitt 
[4] (see also the subsequent works carried in Refs. [5]) in 
a wide variety of situations and with the aid of quite 
different approaches (see Refs. [6] for excellent reviews 
on the subject). Particularly, perturbative and numeri- 
cal approaches were applied for single mirrors [7-10] and 
cavities [11]. Initial field states different from vacuum 
were also considered for single mirrors [12, 13] and cav- 
ities as well [14]. The first experimental observation of 
this phenomenon was recently announced in Ref. [15]. 

Taking into account the difficulties in generating ap- 
preciable mechanical oscillation frequencies (of the order 
of GHz) to obtain a detectable number of photons, recent 
experimental schemes focus on simulating moving bound- 
aries by considering material bodies with time-dependent 
electromagnetic properties. These possibilities were first 
proposed by Yablonovitch [1] and have been further de- 
veloped in theoretical works which considered materials 
with time-dependent permittivities and time-dependent 
surface conductivities [16-18] (see the nice compilation 
done in Ref. [17]). For instance, in Ref. [18] the DCE 
for a massless scalar field within a cavity containing a 
thin semiconducting film with time-dependent conductiv- 
ity and centered at the middle of the cavity was studied. 
The coupling of such a film to the quantum scalar field 
was modeled by a delta-potential with time-dependent 
strength. A generalization to the case of an electromag- 
netic field was carried out in [17]. Very promising and in- 
genious experimental set-ups to simulate non-stationary 
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boundaries include the changing of the reflectivity of a 
semiconductor by the incidence of a periodic sequence 
of short laser pulses [19] or by using a coplanar waveg- 
uide terminated by a superconducting quantum interfer- 
ence device (SQUID). Applying a variable magnetic flux 
on the SQUID, a single moving mirror can be simulated 
[20, 21]. A first step toward the experimental verifica- 
tion of the DCE was recently made in [22] using this 
approach. Moreover, the same group recently claimed to 
have observed the DCE [15]. 

Key ingredients in the predictions of the DCE are the 
boundary conditions (BC) under consideration and nat- 
urally the quantum field submitted to those BC. Quite 
general BC are the so called Robin ones which, for the 
case of a scalar field in 1-1-1 dimensions and a single 
mirror fixed at a: = a, are defined by <j)[t,x — a) — 
^[dx4>it, x)]x=a, where 7 is a real parameter (called here- 
after as Robin parameter). For the case of a moving 
boundary, the previous relation is imposed in the comov- 
ing frame and the corresponding BC in the laboratory 
frame is obtained after an appropriate Lorentz transfor- 
mation. 

This BC have the nice feature of interpolating contin- 
uously Dirichlet (7 — >■ 0) and Neumann (7 — > 00) ones 
and occurs in several areas of physics and mathemat- 
ics. For instance, in classical mechanics they will appear 
if one considers a vibrating string coupled to a spring 
which satisfies Hooke's law and is localized at one of its 
edges [9, 10, 23]. In non-relativistic quantum mechanics, 
Robin BC occur as the most general BC imposed by a 
wall ensuring the hermiticity of the hamiltonian as well 
as a null probability flux through it [24]. Regarding the 
static Casimir effect [25], it was shown that the Casimir 
force between two parallel plates which impose Robin 
BC on a real scalar field may have its sign changed if ap- 
propriate choices are made for the corresponding Robin 
parameters of each mirror [26]. Such kind of repulsive 
Casimir force was also predicted, in the case of paral- 
lel plates, by Boyer in the 70s, who considered a pair of 
perfectly conducting and infinitely permeable plates [27]. 
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Further investigations on the influence of Robin BC in 
the static Casimir effect, including thermal corrections 
and the case of Casimir piston setups, were carried, for 
instance, in Refs. [28]. See also Refs. [29] for the influ- 
ence of this BC on the structure of quantum vacuum. 

Only recently Robin BC were considered in the context 
of the DCE. For a massless scalar field in 1+1 dimensions 
submitted to a Robin BC at a single moving mirror the 
radiation reaction force on the moving mirror and the 
particle creation rate were computed in Refs. [9, 10]. 
Interestingly, for Robin BC, the radiation reaction force 
acquires a dispersive component, in sharp contrast with 
Dirichlet and Neumann cases where the force is purely 
dissipative. It was also shown that, for a given Robin pa- 
rameter, there exists a mechanical frequency of motion 
that dramatically reduces the particle creation effect [10]. 
Finally, and of crucial importance for the present work, 
Robin BC can also be useful to describe phenomenologi- 
cal models for penetrable surfaces and under certain con- 
ditions they simulate the plasma model for real metals 
[30]. In these situations, for frequencies w much smaller 
than the plasma frequency wp, the Robin parameter 7 
can be identified as the plasma wavelength Ap. In other 
words, the Robin parameter 7 gives us an estimative of 
the penetration length of the mirror under consideration 
[32]. 

Since to simulate a motion of a reflecting mirror is 
equivalent to simulate a real metal with time-dependent 
plasma wavelength, the above interpretation of 7 leads 
naturally to the consideration of time-dependent Robin 
parameters. Specifically speaking, it is quite natural to 
simulate the motion of a reflecting mirror by considering 
the quantum field submitted to a Robin BC at a static 
mirror but with a time-dependent Robin parameter "f{t). 
The kind of boundary motion which is being simulated 
is determined by the kind of time-dependence of jit). 
The purpose of this paper is precisely to analyze this 
situation for a massless scalar field in 1+1 dimensions. 
Particularly, we shall compute explicitly the particle cre- 
ation rate for a natural choice of time-dependence for 
7(t) which is directly related to recent experimental pro- 
posals. This paper in organized as follows: in Sec. II the 
Bogoliubov transformation between the in and out cre- 
ation/annihilation operators are obtained, allowing us to 
find the spectral distribution of the created particles and 
the particle creation rate in Sec. Ill and IV, respectively. 
Finally, in Sec. V we present our conclusions and final 
remarks. Throughout this work we consider h = c = 1. 



II. THE BOGOLIUBOV TRANSFORMATION 



tion, d'^(j> = 0, and is submitted to a time-dependent 
Robin BC at a mirror fixed at the origin, namely, 
'j{t)d(j)/dx\x=o — ^(0,i) = 0. For simplicity, we assume 
that 7(t) departs only slightly from a positive constant 
7o, so that we can write ^{t) = 70 + 5j{t), where 5^(1) is 
a smooth time-dependent function satisfying the condi- 
tion max|(57(t)| ^ 70, for every t. Under these assump- 
tions in the limit 70 — > c» we recover Neumann BC. On 
the other hand, to reobtain Dirichlet BC (70 0), be- 
cause of condition max | (57(f) | <gC 70, we must also take 
5'y{t) — > 0. If we consider only 5'y{t) = we re-obtain the 
usual time-independent Robin BC. Moreover, we shall 
also impose that 6j{t) foY t ^ ±00. The BC satis- 
fied by 5'y{t) then reads 



To 



d(j){x, t) 



dx 



- 0(O,O+^7(i) 



x=0 



d(j){x, t) 



dx 



x=0 



(1) 



Also for the field, a perturbative approach will be 
adopted. Following Ford and Vilenkin [31] we write 



(t>{x, t) = (j)o{x, t) + 54>{x, t) , 



(2) 



where, by assumption, (j)o satisfies the Klein-Gordon 
equation, d'^(po = 0, and the time-independent Robin BC, 



7o 



d(f>o{x,t) 
dx 



(/>o(0,t) = 0. 



(3) 



x=0 



The small perturbation 5(j) takes into account the contri- 
bution to the total field (f) caused by the time-dependence 

of the Robin parameter, described by the function S"f{t). 
Since both (p and 4>o satisfy the Klein-Gordon equation, 
so does 5(f), namely, d'^5(f) = 0. The BC satisfied by Sep 
is obtained, up to first order terms, by substituting (2) 
into Eq. (1), which leads to 



70 



d5(t){x, t) 



dx 



6cj>{0,t) = -6^{t) 



x=0 



d(j)o{x,t) 



dx 



x=0 



(4) 



We start considering a real massless scalar field (p in 
1+1 dimensions which satisfies the Klein-Gordon equa- 



where Eq. (3) was used. Hereafter it will be convenient 
to work in the Fourier domain, such that 
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^{x,Lj) = j dt (l){x,t) e'"^* ; ^o{x,uj) = j dt (t)o{x,t) 



iijjt 



(5) 



r 



It is worth emphasizing at this moment that, by assump- 
tion, (57 is a prescribed fimction of t, so that 5V{uj) is 
known, in principle. Since (j)o{x,t) is the solution with 
time-independent Robin BC, this field is already known, 
and so docs its Fourier transform, which is given by (for 
the region a; > 0), 



^o{x,uj) 



4tt 



|a;|(l+7o^c.2) 



sin(wa;) + 70C1J cos(wa;) 



X e{uj)a{u>) - Q{-u>)a\-uj) 



(6) 



where 0(a;) is the Heaviside step function. The operators 
a{uj) and a^(ix)) satisfy the usual bosonic commutation 
relation [a(w), a^(w')] = 2ttS{lo —oj'). 

In order to obtain ^{x.w) = ^q{x,uj) + 5^{x,lj) we 
need to compute 5^{x,lj), which satisfies the Helmholtz 
equation. 



(dl +uj^) 5^{x,uj)=Q, 



(7) 



and is submitted to the BC below, obtained by Fourier 
transforming Eq. (4), 



70 



f doj' 

J 



d6^x, 

dx 
d<^o{x,uj') 
dx 



(5$(0,w) 



x=0 



(8) 



x=0 



A further condition that must be imposed to the solution 
of Eq. (7) for X > is that it will lead to a solution 
for (j){x,t) that must travel to the right, since 5(j){x,t) 
must describe a contribution coming from the mirror, 
and not going towards the mirror. The desired solution 
can be written in terms of Green functions. Following 
the procedure given in [10] it can be shown that the in 
and out fields, denoted respectively as $in and ^out, are 
related to each other according to 



$out(a;,a;) = $in(a;,w) + — \G'^^{Q,x,uj) 

7o L 



-Ggi-(0,a;,u;) 



X i7o 



d5^{x,u:) 



dx 



5$(0,w) 



x=0 



(9) 



where G'^\Q,x,uj) (G^^^(0,a;,a;)) is the retarded (ad- 
vanced) Robin Green function, satisfying the time- 
independent Robin BC at x = 0. These Green functions 



are given, respectively, by 



and 



G^'l-(0,x,a;) 



70 



1 — i'yooj 



70 



■ i7ow 



(10) 



(11) 



Inserting Eqs. (6) (appropriately relabeled as <I>out and 
(f>in), (8), (10) and (11) into Eq. (9), we can readily 
obtain the Bogoliubov transformation between aout and 
ttin and its hermitean conjugates: 



aout(w) = ain(a;) - 2i,l ^ 

V 1 + 70^^ 



+ 00 



duj' 



v2,,,/2 



) 27r y 1 -H 75a;' 
e(w')ain(w') - @{-u:')aU-uj')\ 5V{uj - J). 



(12) 



Noting that the annihilation operator aout(i^) is given in 
terms of the annihilation and creation operators a\^{Lu) 
and a\^{u)), respectively, we conclude that the state |Oin) 
is not annihilated by the aout('^) operators. Conse- 
quently, we can state that particles were created from 
an initial vacuum state due only to the time-dependence 
of 5^{t) in the BC (1) imposed on the field by the 
static mirror. In fact, for 5^{t) — for all times, 
which corresponds to a static mirror imposing the stan- 
dard time-independent Robin BC on the field, we have 
flout (^) = O'lni^) and no particles will be created, as 
expected. The particle creation effect will be further in- 
vestigated in the next sections, where we will choose a 
specific time-dependent expression for 7(t) in order to 
compute explicitly the corresponding spectral distribu- 
tion of the created particles as well as the respective par- 
ticle creation rate. 



III. SPECTRAL DISTRIBUTION OF THE 
CREATED PARTICLES 

We start by writing the spectral distribution of the 

created particles as 



dN{uj) 



duj 



1 



(OinI Clout ('^)«out(w) |Oin>rft^, (13) 



duj " 2tt 

where dN{ui)/duj is the number of created particles with 

frequency between lu and w + dcj (cj > 0) per unit fre- 
quency. From the previous definition for dN{w)/du}, it 
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follows immediately that the total number of created par- 
ticles from t — —oo to t = +00 is given by 



N 



dN{uj) 
duj 



duj . 



(14) 



From Eq. (12) and its hermitian conjugate 0^^^., 
straightforward to show that 



it is 



dN{uj) 
duj 



-72^2 



duj' 



27r 1 



x\5V{l^J - l^J' 

)l e(^')- (15) 



In what follows we will obtain the spectral distribution 
for a particular case of 5T{uj). With this purpose in mind, 
let us consider the following expression for 5^{t)^ 



<57(i) =eocos(^^o0e"l*l/^, 



(16) 



with loqT 3> 1. This choice of 5^{t) may simulate, for 
instance, the changing magnetic flux through a SQUID 
fixed at the extreme of a unidimensional transmission 
line, as in Ref. [20], where a Robin-like BC arises natu- 
rally from quantum network theory applied to the system 
under consideration. 

The expression of 5T{uj), obtained by Fourier trans- 
forming Eq. (16), contains, in the limit of uj^T ^ 1, two 
sharped peaks around w = iwo, which can be approxi- 
mated by Dirac delta functions, leading to the result 

\5V{Lo f « ^elT[5{u: - c^o) + 5[u + u^)] . (17) 

Substituting the above result into Eq. (15), we finally 
obtain the desired spectral distribuition. 



dN(L0) 

duj 



2tt 



UJ {lUq — Uj) 



(l+7o'^2)[l+7o'(^o-^)^ 



e(c^o 



(18) 



for this particular situation. 

A few comments are in order. Firstly, observe (see Fig. 
1 and Eq. (18)) that dN{u))/du) vanishes for w > wq, 
which means that no particles are created with frequen- 
cies larger than wq — the characteristic frequency of the 
time-dependent BC. We also notice that the spectrum is 
left invariant under the replacement w — > wq — w. This is 
a signature of the fact that particles are created in pairs: 
for each particle created with frequency uj there is a twin 
particle created with frequency wq — In second place, 
note that for eq 0, where a Robin BC with a time- 
independent parameter 70 is re-obtained, the spectrum 
of created particles vanishes, as expected (recall that the 
mirror which imposes the BC on the field is at rest). Fur- 
ther, for a fixed (finite) value of uq, the limit 70 — > 00 
(Neumann BC imposed on the field at a static mirror) 
also leads to a vanishing spectrum of created particles. 
Finally, since we assumed eq ^ 70i the limit 70 — > 
(Dirichlet BC imposed on the field by a static mirror) 
necessarily leads to a vanishing spectrum as well. 



i5 0.1 




0.05 



FIG. 1: (Color online) The spectral distribution of the created 
particles [(27r) / (eq T)\ dN/duj as a function of cj/ojo for several 
values of 70. Notice the reflexion symmetry around cj/cjq = 
0.5: a signature of the fact that particles are created in pairs. 
The full line corresponds to 70 = 1; the dashed line to 70 = 5 
and 20 x [(27r)/(eo T)] dN/dw; and the dotted line to 70 = 10 
and 100 x [(27r)/(eg T)] dN/duj. 



IV. PARTICLE CREATION RATE 

The total number of created particles is obtained by 
substituting Eq. (18) in (14), namely. 



N = 



2tt 



Ztt 



UJ [lJq — LO) 0(wo ^ ^) 

(l+7o^')[l + 7o(^o-w)2 



dhJ 



where f = 7otj-'o and the function i^(^) is given by 
(2 + e^) In (1 + ^2) _ 2^ arctan(0 



(19) 



(20) 



As N is proportional to T - as expected for an open cavity 
- the physical meaningful quantity is the particle creation 
rate defined as i? = N/T, that is 



R 



2 3 
2tt 



(21) 



In the limits 70^9 ^ 1 and 701^0 3> 1, the particle cre- 
ation rate are approximately given by 



4^0 

127r 



for 70^0 < 1 (22) 



2n 



21n(e) 



for 7oa;o > 1. (23) 



For the sake of comparison with Eq. (21), we recall 
the total particle creation rates for moving mirrors with 
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static mirror witti time-dependent Robin BC 
I\/loving mirror witti Robin BC 
Moving mirror with Dirichilet/Neumann BC 
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FIG. 2: (Color online) Comparison between the total particle 
creation rates given by Eqs. (21), (24) and (25). The full line 
corresponds to scaled creation rate 10 x [(27r7o) / (eo)] ^• 
The dashed line corresponds to 10 x [(27r7o) / ((^go)] ^ti-R- 
Finally, the dotted line corresponds to [(27r) / ((Jqo)] ^d/N- 
In both curves involving Robin BC we considered 70 = 1. 



FIG. 3: (Color online) Behavior of Eq. (21) for a larger 
range of ujq, assuming 70 = 1. The creation rate approaches 
smoothly to zero as ijjq ^ 00. This does not happen for the 
other cases showed in the Fig. 2: for movmg mirrors the par- 
ticle creation rate increases unlimitedly for larger values of 

LOO- 



Dirichlet [8] (or equivalently for Neumann BC as proved 
in [13]) 



Rt 



and for time-independent Robin BC [10] 



-Rti-R 



^ ) G(7o^o), 



(24) 



(25) 



where 



S. [4^ + + 12 arctan(^)] - 6 (2 + In (l + 



ee (4 + e 



(26) 



The formulas above were obtained assuming a non- 
relativistically small amplitude oscillatory law of motion 
for the mirror. For both cases Sqq is the amplitude and 
ujQ is the frequency of oscillation. We remark that for 
jqujq <^ 1 the particle creation rate in our model is ex- 
actly the same as for of a moving mirror [8] with Dirichlet 
BC where eq plays the role of the amplitude of oscilla- 
tion of the motion. This reinforces the possibility of sim- 
ulating moving boundaries through a static mirror with 
time-dependent Robin BC. The three particle creation 
rates are compared in Fig. 2. 

It is worth noting that the particle creation rate shown 
in Fig. 3 starts growing with ujq until it achieves a maxi- 
mum value for a given value of and then it approaches 
monotonically to zero as loq goes to infinity. This be- 
haviour should be compared with that obtained for a 
moving mirror which imposes on the field a Robin BC 
with a time-independent parameter, where the particle 
creation rate after passing through one maximum and 



one minimum grows indefinitely as goes to infinity 
(see Ref. [10]). Naively, we could expect similar behav- 
iors for these two problems, after all, a time-dependent 
Robin parameter should simulate, in principle, a moving 
mirror so that a high frequency oscillating ^(t) should 
mean a high frequency oscillating mirror. However, the 
interpretation of the Robin parameter 7 as an estimative 
of the penetration depth of the material boundary is rig- 
orously proved only for static mirrors. Even in this case, 
this identification is valid only for the field modes whose 
frequencies are much smaller than the plasma frequency 
(but this condition is easily achieved since the plasma 
frequency is much higher than the mechanical frequen- 
cies we want to simulate). It is plausible that such an 
interpretation remains valid for slowly time- varying j{t), 
but not for high frequency oscillating ^{t). In fact, our 
results show that this interpretation for ^(t) fails for high 
values of ujq. 



V. CONCLUSIONS AND FINAL REMARKS 

Exploring the peculiar properties of Robin BC, par- 
ticularly, the interpretation of the Robin parameter, we 
presented a simple and yet instructive theoretical model 
where a single static mirror with time-dependent prop- 
erties described by a time-dependent Robin parameter 
simulates a moving boundary. We used this model to 
study analytically the dynamical Casimir effect of a sys- 
tem that may of some value for further understanding of 
a ongoing experiment based on a one-dimensional trans- 
mission line terminated by a SQUID. In this setup a time- 
dependent magnetic flux through the SQUID gives rise to 
particle creation phenomenon. Employing a perturbative 
approach, we showed that particles can be created due 
to the time-dependence of the Robin parameter 7. We 
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obtained explicitly the spectrum of the created particles 
as well as the total particle creation rate for a particular 
choice of jit) which has a practical interest concerning 
the experiment just described. Our model can also be 
used as a theoretical model to investigate other exper- 
imental setups suggested for measuring the dynamical 
Casimir effect, as for example, the promising experimen- 
tal proposal of the Padua group [19]. All we have to do 
is to choose appropriately the time-dependence of 7(i) to 
simulate correctly the physical situation under consider- 
ation. 

We emphasize that the particle creation phenomena 
due to a time-dependent Robin BC imposed on the field 
at a static mirror has similarities and differences with the 
case where a time-independent Robin BC is imposed on 
the field at a moving mirror, as discussed by Mintz et al 
[10]. The main difference being the respective behaviors 
of the total particle creation rate for high values of ojq 
(in the former case, where ojq means the mechanical fre- 
quency of the moving mirror, this rate grows indefinitely 
as Wo — >■ 00, while in the latter case, where wq gives a 



measure of how quick the time-dependent Robin param- 
eter "f{t) varies, this rate goes to zero, as loq — > oo). In the 
appropriate limits of the usual time-independent Dirich- 
let (70 — >■ 0), Neumann (70 00) and Robin {Sj{t) = 0) 
BCs no particles are created, as expected. The general- 
ization of the present work for 3-1-1 dimensions and cav- 
ities are also expected to have induced photon creation. 
These issues are under investigation and will be discussed 
elsewhere. 
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